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Abstract 

The spin-extended semiclassical chiral fermion (we call the S-model), which had been used to 
derive the twisted Lorentz symmetry of the “spin-enslaved” chiral fermion (we call the c-model) 
is equivalent to the latter in the free case, however coupling to an external electromagnetic 
field yields nonequivalent systems. The difference is highlighted by the inconsistency of spin 
enslavement within the spin-extended framework. The S-model exhibits nevertheless similar 
though slightly different anomalous properties as the usual c-model does. The natural Poincare 
symmetry of the free model remains unbroken if the Pfafhan invariant vanishes, i.e., when the 
electric and magnetic fields are orthogonal, E ■ B = 0 as in the Hall effect. 
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1. INTRODUCTION 


The semiclassical chiral model (we call here the c-model) allows for a derivation of the 
chiral magnetic effect and the chiral anomaly, respectively, bypassing complicated quantum 
calculations HHg. The free c-model, which has no genuine spin degree of freedom, carries a 
curious “twisted” Lorentz symmetry [6Hin], conveniently derived by relating it to Souriau’s 
massless spinning particle m- The latter (we call the S-model), carries a mass-zero, spin-s 
Poincare symmetry. Compared to the c-model, the S-model has two additional degrees of 
freedom represented by an “unchained” spin vector, s, whose projection onto the momentum 
is fixed, s ■ p = s PE]. Free spin can however be “enslaved” to the momentum, 

s = sp, (1.1) 


by a suitable “Wigner-Souriau translation”, which embeds the free c-model into that of 
Souriau m, making them equivalent PEj. 

The c- and S-models are no longer equivalent, though, when the systems are put into a 
field, as highlighted by the explicit solution presented in Sec. 5 A of p. In particular, spin 
can no longer be consistently enslaved within the S-model p. 

In this Letter we show that the minimally coupled S-model, although nonequivalent to the 


c-model, admits nevertheless similar transport properties, see eqns ( |5.3[ )-(6.3)-(6.5) below. 

We formulate our results within the framework of Souriau pm. To make our paper more 
self contained, we remind the reader of some basic facts while referring to these references 
for details. 

The ultimate description of a mechanical system is provided by its space of motions, M, 
which is a symplectic [and therefore an even-dimensional] manifold; its symplectic form, 
u, is regular. A symmetry group acts on M by preserving its symplectic form, uj. If the 
symplectic action is transitive, then (M, ca) is identified with a coadjoint orbit carrying its 
canonical symplectic form. 

Conversely, one can start by constructing such an orbit and then seek a physical inter¬ 
pretation for it. For this end, it is convenient to use what Souriau calls an evolution space, 
V, which is endowed with a closed two-form a of constant rank r > 1. The distribution 
provided by ker a is integrable, and its characteristic leaves are identified with the classical 
motions of the system. The space of motions is the quotient of the evolution space by the 
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characteristic foliation of cr, 


(M, a;) = (y, (j)/ker cr. (1.2) 

Thus dim V = dim M + r. The points of M are be labeled by constants of the motions. 

The concrete realization of this abstract framework in our context is summarized in 
section |3] below. 


2. THE C-MODEL 

We hrst summarize some aspects of the c-model we will need to be compared with 
those in the S-model. In Souriau’s framework [siiiiin] the model can be described by 
a 7—dimensional evolution space with coordinates {x,p,t), endowed with the closed 
2—form ac = oJc — dh A dt with one-dimensional kernel, where the symplectic form and the 
Hamiltonian are, 


ujc = dpi A dxi + -CijkBidxj A dxk 


1 


A\p\' 


t^ijkPidpj A dpfc. 


h= \p\+' 


respectively, where — V0 = E [HKH]. The Poisson brackets are therefore 


{Xij Xj^ ^ijk~- 


{Xi,Pj}= 


6{j -j- cBjhj 


{PiiPj} ^ijk: 


eBj, 


( 2 . 1 ) 


( 2 . 2 ) 


1 + eb B' ' 1 + ebB' ^ 1 + eb ■ B 

where b = sp/|pp, p = p/|p| is the “Berry monopole” of strength s in the momentum 
space [12]. Usually s = 1/2 [THS|. In the denominator we recognize here the square-root of 
the determinant of the symplectic matrix Uc, 1 + eb ■ B = y^det(ciyj = Dc- The system is 
regular when Dc ^ 0. The Hamilton equations 


dx 

Dc — = p + eE X b + e{b ■ p)B, 
Dc^ = eE + ep X B + e^{E ■ B)b, 


(2.3a) 

(2.3b) 


reproduce eqns. (14)-(15) in [T] and correspond to the one-dimensional kernel of dc j6l fTn[T2] . 
Factoring out the kernel yields the 6 dimensional space of motions (identihed here with the 
phase space) [HI [H]. Particular solutions putting in evidence the role of the anomalous 


velocity in (2.3) were studied in 


The invariant volume element is the 6/2 = 3rd power of the symplectic form Uc [HI [I2| 
and its pull-back to the evolution space is the 3rd power of cxc. 


dVc = Dcd^pd^x. 


(2.4) 
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Then Liouville’s theorem takes the anomalous form 


dDr d(Dcx) d(Dcp) 

+ r. + =(E ■ B) Vr 


p 


2|p| 


= 2we^(E- B)6'Hp). 


(2.5) 


dt dx dp 

Let f{x,p, t) be a distribution on the phase space which we assume to satisfy the collision- 
less Boltzmann equation dtf + dxf x + dpf p = 0. The current density is, 


j = / fxDc 


d?p 

(27r)3 


= fP 


d?p 

(27r)3 


+ eB 


f d^P 

2|p|2 (27r)3 


+ eE X 


fP d^P 
2|pP (27r)3 


• ( 2 . 6 ) 


The hrst term on the r.h.s. is the normal current, the second one represents the chiral 

magnetic effect (CME) and the last one is the anomalous Hall current [T|. Dehning the 

f d^p 

particle density as p{x^t) = / yields the anomalous continuity equation (referred 


(2^) 


to as the chiral anomaly), 
where /o is the value of the distribution function at p = 0 m- 


dtp + V-j = ^{E-B)fo, 


(2.7) 


3. THE MASSLESS SPINNING MODEL, MINIMALLY COUPLED TO AN E.M. 
FIELD 

The evolution space of the S-model, U®, has, w.r.t. to the c-model, two additional degrees 
of freedom represented by the spin vector, s, whose projection onto the momentum is fixed, 
s ■ p = s [6]. However, the kernel of the free two-form cxc which yields the equations of 
motion is now ^-dimensional spanned by Wigner-Souriau translations, yielding, once again, 
a 6-dimensional space of free motions iniiii]. 

It is worth, at this point, to compare the c- and S-models. In the c-model the evolu¬ 
tion space, U^, is 7-dimensional and the motions are curves, which are tangent to the 1- 
dimensional kernel of cXc. In the S-model instead, the evolution space, V^, is 9-dimensional; 
the kernel is as is 3 dimensional. The respective spaces of free motion are, therefore, 6 
dimensional in both cases. The key point is that factoring out the respective motions yields, 
in both cases, the same space of motions, (M,cn), as illustrated on Fig. 

Coupling to an external electromagnetic field is introduced through Souriau’s minimal 
coupling schema which requires to add to the free form a [e-times] the electromagnetic 
held tensor, written in terms of the “true” position, r, - the one which transforms in the 
usual way under a Lorentz boost [6]. 
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FIG. 1: The free motions of the c-model are described by curves in the 7-dimensional evolution 
space, whereas the motions of the S-model are 3-dimensional surfaces lying in 9-dimensional evo¬ 
lution space G®. Factoring out the motions yields, however, the same space of motions, 
for both systems. 

The S-model is thus described by a 9-dimensional evolution space with coordinates 
r,p7^0,s,(s-p = s = i), endowed with the closed 2-form as [6]. The Hamiltonian is still 
of the form h = \p\ -\- ecf, remember however that the potential (assumed static) is now a 
function of the “true” position r, (f = 4>{r). Spelling out in 3 -|- 1 dimensional notations the 
constraints which define the spin-extended evolution space and the symplectic form, 

= 0 , = 0 , 


(3.1a) 

^I\ pU.Lj ^-r (3.1b) 

cf. eqns. # (5.3) and # (3.4) in [U], a long but straightforward calculation using also the 


a = -dP^ A ^ dS\ A dS^^ + dR^^ A dF^ 


constraints in eqn. (3.1a) yields the complicated-looking symplectic form as = ujs ~ dhAdt, 


uis = dpi A dvi P -CijkBidrj A dr*, 


IpI^ 


[eijkSi + 2pj{p X s)k]dpj A dpk 


Pi 

-h 2[eijk{^ - Si) + 2{p X s)jPk]dsj A dsk 
2 1 

- Tl [(^ijkis'^Pi + vSj) + (p X s)jSk + {pj - Sj){p X s)k]dpj A dsfc. 


(3.2) 


IPI — ' - 2 

If spin could be enslaved, s = sp as in the free case, terms with p x s would drop out. As 
we show below however, spin enslavement is consistent with the S-dynamics only in the free 


case in that (5.1) would reduce to (2.1) in the free, but not in the coupled case, when the 


two models are radically different. 
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The rather weirdly-looking equations of motion calculated from the kernel of as, 


dv 

{p-B)- = B-pxE, 

(P- = e{E ■ B)p, 

ds 

(P'-B)— = px B - px {px E), 


(3.3a) 

(3.3b) 

(3.3c) 


are similar to but different from the analogous equations, (2.3), in the c-model. Note in 


particular the absence of the usual momentum on the rhs. of the velocity relation which is, 
so to say, “purely anomalous”. The upper two equations here are decoupled from the lowest 
one, so that the space-time motion does not depend on the spin at all. We record for later 


use that the direction of the momentum is a constant of the motion, dp/dt = 0. Eqn (3.3c) 


shows, moreover, that s = p x r. Thus, although spin is not more enslaved, its motion is 
entirely determined by that in space-time. 


The eqns (3.3a-c) are valid under the regularity assumptions [U] 


{i) p ■ B ^ 0, (a) s ■ {B — p X E) 0. 


(3.4) 


These conditions are preserved by the dynamics, as shown by using the equations of motion. 
The hrst of these conditions will be interpreted below as the non-vanishing of the system’s 
determinant. 

It is worth stressing that the neither the free equations of motion nor the free Souriau 
form can be recovered by simply letting the fields go to zero. This limit is in fact a 


singular one as seen from eqn. (3.4) and is highlighted by the fact that turning off the fields 


converts the one-dimensional motion curves into 3-dimensional surfaces, cf. [B]. 

As proved in [6] and seen also directly, the helicity condition s ■ p = | is consistent with 
the coupled equations of motion. 


4. “DYNAMICAL” POINCARE SYMMETRY 

Before turning to study the transport properties, we would like to point out a rather 
curious fact. For the free S-model, the angular momentum vector, 

£ = r X p + s, (4.1) 
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is plainly conserved [ 6 ]. Now for arbitrary constant electric and magnetic fields, E and B, 
the eqns of motion imply that various terms cancel, leaving us with 


d£ (E-B) 

-n = e ^ ^ XP, 

at p ■ B 


(4.2) 


[assuming the determinant does not vanish, p ■ B 7 ^ 0]. Therefore (4.1) is not conserved in 
general, as expected. A surprising observation is, though, that when E and B are orthogonal 
so that the Pfaffian invariant vanishes, —\{i<F.F) = —E ■ B = t), then all three components 
of the angular momentum are conserved and the system has a full rotational symmetry. 

For comparison, the conserved angular momentum in the c-model is = '^ x p + i.e.. 


the spin contribution is enslaved to the momentum. Our calculation leading to (4.1) shows. 


however, that £c is not conserved in the minimally coupled S-model, £c 7 ^ 0, and it is the 
“unchained component” s — | p of spin which is restores angular momentum conservation. 
Further aspects of the angular momentum for are reviewed in |14j . 

The unbroken rotational symmetry in Hall-type crossed e.m. helds we found here above 
can actually be extended into a full Poincare symmetry [12]. The equations of motion 


(3.3a c) imply, in 4D notations, that the quantity 


+ eFfR^ 


(4.3) 


reminiscent of “magnetic translations” in the massive Landau problem is conserved, and a 
short calculation shows that 


p/^ = and , 


(4.4) 


2S.F 2S.F 

where ^ is the generalized Pauli-Lubanski vector [T3|. Therefore when 

the Pfaffian vanishes, both the linear and the Lorentz momenta are constants of the motion. 


P = const and M = const if -kF.F = 0 (4.5) 

i.e., the full Poincare momentum (M, P) is conserved extending what we had found for the 
angular momentum. 


5. TRANSPORT PROPERTIES 

Turning to the transport properties, we choose .B 7 ^ 0 to point into the 3rd direction, 
B = B z and eliminate one component of the spin vector s, say -Sa = ^(| — Sipi — S 2 P 2 ), 
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leaving us with the 8 independent coordinates r,p, si, S 2 - Then a tedious calculation yields 


Us = dpi A dvi + leijkBidvj A dvk + 


S3 


2\p\p: 


-eijkPidpj A dpk 


— \p 1 p 2 dp 1 A dsi + (1 — Pi)dpi A ds 2 — (1 — p 2 )dp 2 A dsi] 

P3 

— \ - Plp 2 dp 2 A dS 2 + p2p3dp3 A dsi - Pip 3 dp 3 A dS 2 ] . 

P3 


(5.1) 


Note that no dsj A dsk terms show up in (5.1). Then a lengthy calculation yields the 
determinant of the symplectic form, 

ep■B eB 


Ds = \/det(ws) = 


\P\‘^P3 


\P\ 


(5.2) 


This result is surprisingly simple and somewhat unexpected in that it does not involve the 


spin. Note that, compared to (2.4), (5.2) has a “naked” b ■ B term, but the normal “1” is 
missing. 

When 7 ^ 0 the system is regular; the motions follow curves (world lines) so that the 
space of motions is 8 -dimensional with coordinates (r,p, si,S 2 ). 

What happens when p ■ B —)■ 0 ? The determinant goes to zero, Ds —)■ 0, and the 
system becomes singular, necessitating reduction, analogous to what happens for an “exotic 
particle” in the plane |T 6 ]. The characteristic world lines degenerate into three dimensional 
world sheets : the space of motions of the free massless spinning particle is 6 dimensional, 
with coordinates (r,p) alone [HI 171 ITT]. 

In the regular case Ds 7 ^ 0 we assume henceforth, the Liouville theorem takes now the 
form (a = 1, 2) [15], 


dDs , d{Dsr) , d{Dsp) , d{DsSa) 


dt 


+ 


dr 


+ 


dp 


+ 




dSa 


P\^P3' 


^2 (.E^ B) 

P3 


This result only differs from (2.5) in a factor 2 and in the factor p^^ which is in fact a 
constant of the motion, as we noted earlier. Paradoxically, the right hand side here, in this 
spin-extend model, is independent of the spin, as long as the latter does not vanish - whereas 


the “spin-enslaved” c-model (2.5) has a spin-remnant, namely the factor s = 1/2. [20] . 

It can be inferred from the equations of motion that the helicity constraint p ■ s = s 
is preserved by the S-dynamics. However the explicit solutions found in [ 6 ] indicate that 


enslavement, s = sp in ( 1 . 1 ), is not preserved : spin ean not be eonsistently enslaved within 
the minimal S-model which is therefore dehnitely different from the c-model. This is also 


















obvious by counting the degrees of freedom: the electromagnetic held breaks the Wigner- 
Souriau translations and reduces the dimension of the kernel from 3 to 1, therefore the space 
of motions is 8, and not 6-dimensional. Unlike as in the free case, spin is a genuine degree 
of freedom, which can not he eliminated. 

In the non-singular case yet another tedious calculation allows us to hnd the Poisson 
brackets, 


{ri,rj} 

{si,rj} 

{S1,S2} 


^ijkPk 
ef) ’ B’’ 




{ShPj} 


S3 - 


|pb3 

ep ■ B 


Bjpj 

p-B' 


{pi,Pj} = 0 , 


\P\ 

p ■ B 


{eijkBk +Pi(p X B),), 


(5.4) 


which are substantially different from those for the c-model, (2.2). The Jacobi identities 
follow from dus = 0, and can also be checked directly. Note here the absence of the usual 
“Heisenberg” term 6ij in {ri,Pj}, similar to the dropping out of the momentum term p from 


the velocity relation in (3.3a). Note also that the momenta commute instead of closing 
on the magnetic held, as expected. It is therefore reassuring that the associated Hamilton 


equations yield (3.3a)-(3.3b)-(3.3b) as they should. Thus, the diherence between the coupled 
c- and S-models originates in both the symplectic structure and the Hamiltonian. 


6. CHIRAL MAGNETIC EFFECT AND CHIRAL ANOMALY 


The particle current is determined in terms of the coordinates r,p,si,S 2 using the de¬ 


terminant of the symplectic matrix (5.2). The invariant phase space volume element of the 


8-dimensional space of motions is Vs = = a\/A\ [HI [12] be., by (5.1), 

eB 

dVs = Ds dPrdPpdsids 2 = i —— dPr dPpdsids 2 

\P\ 

also expressed in a covariant way useful for calculations. 


(6.1) 


dVs = ^ ■ s — dPrdPpdPs . 

IpP ^ 2^ 


( 6 . 2 ) 


If f{r,p,si,S 2 ) is a distribution function on the spin-extended space of motions, which 
we assume to satisfy again the collisionless Boltzmann equation, now —l-ri 7 r- \-piT: -h 

Ot OUi opi 
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df 

SaT :— = 0 (* = 1,2,3, a = 1,2), the the particle current, 

OSa 


fP 


j{rA) = [fr Dsd^p dsids 2 = eB [-— l—d^pdsids 2 + eEx [ , d^pdsids 2 (6.3) 

J J \prP3 J \prP3 

is decomposed into merely two (and not three) terms, namely into a chiral- magnetic and an 
anomalous current analogous to those, ( 2.6[ ), for the c-model [Il|2]. The absence of a normal 
current follows from that of the usual p. The particle density p{r,t) = J f Dsd^p dsids 2 , 
satisfies 

dp{r,t) f df 


dt 


d D P d J' 

(^-^'jf(fpdsids2+ / Ds-fd^pdsidsi- 


(6.4) 


Dropping boundary terms we hnd, for constant helds B, E and no explicit time dependence, 
^^^ + Vr-j{r,t) = e^E-B [fVp ■ {-^)d^pdsids2 = e^E ■ B (6.5) 

dt J \p\P3 P3 

which differs from ( 2.7[ ) valid for the c-model by the same factors as (5.3) does from (2.5) 

and/o = f{r,p = 0 ,si,s 2 )dsids 2 . 


CONCLUSION 


In this Letter we demonstrated that the spin-extended chiral model (our S-model), in¬ 
strumental in deriving the twisted Lorentz symmetry of the c-model is eqnivalent to the 
latter only in the free case but not when coupling to an external held is considered. This is 
highlighted by the inconsistency with the S-dynamics with “enslaving”. 

The difference comes from the different choice of what one considers as “position”: the c- 


model is coupled to the e.m. held by viewing x in (2.1) as a position [IH3], with no attention 
paid at its “twisted” behavior under a Lorentz boost [61 IBl-fTU]. In the S-model instead, the 
coupling is introduced in terms of the “true” position, r, which does transform in the usual 
way under a Lorentz boost |6]. Let us emphasize that the S-model and its coupling to an 
external held follow from First Principles - namely of Souriau’s Mechanics m- We stress 
that our “true position”, r, is dehned on the evolution space and not on M\ it is the 
combination of position, time, momentum and spin, 

ri V 

(7.1) 
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X = r = r — pt -\ -1—— 

\p\ \p\ 


which is conserved and can label a point of M - i.e., a motion. The Poincare group acts 
naturally on U®, whereas its action projected onto the space of motions is “twisted” and is 
not natural [6]. 
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We note also the the expression (4.1) of the angular momentum “lives” on the evolution 


space. Using the space of motion coordinates x and p allows us to absorb the evolution 
space coordinates into conserved quantities and convert the angular momentum it into 


i = X X p + ^p. 


(7.2) 


In particular the “unchained” part of the spin vector s has been absorbed into x, leaving 
us with the “enslaved” contribution |p. 

The X used in the c-model is in turn a label of c-motions obtained by putting f = 0 into 


x{t) = x+pt, obtained by integrating the c-equations of motion, (2.3), in the free case with 
initial condition x{t) = x. Therefore viewing x in the c-model as a “position,” is, in our 
opinion, unjustihed. 

Returning to the coupled case, we mention that it has been suggested [21 El E] to modify 
the c-Hamiltonian by adding a term. 


h —)■ IpI 


p ■ B 


(7.3) 


Such a modihcation is certainly possible and can be generalized to the spin-dependent case 
[6]. Anomalous coupling yields in fact the dispersion relation 

P 


S = ^\p\^-{eg/2)S-F, S ■ F = s ^ y, E) . (7.4) 

where the real number g represents the gyromagnetic ratio [6] . For s = ^ g = 2 and a weak 


purely magnetic held, (7.4) approximately reduces to (7.3). In this Letter, we studied the 


minimal case g = 0. Our study will be extended to anomalous coupling elsewhere. 

The S-model exhibits properties which are similar to those the chiral one, namely the 
chiral anomaly, CME and AHE. Its advantage is its manifest Lorentz invariance. 

Remarkably, the anomaly vanishes precisely when the system carries a full Poincare 
symmetry - namely when the Pfafhan invariant vanishes. 

We just mention that the additional spin degree of freedom would allow us also consider 


the spin current dehned, by analogy to (6.3), as 


3s = j fsDsd^pdsids 2 

= —e [ , {p ■ E) d^pdsids 2 — eB x f df‘pdsids 2 + eE f d^pdsids 2 - (7.5) 


IPI Ps 


IPI P3 


IPrPs 
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We have here again three terms as in (2.6). Note also the the magnetic and anomalous Hall 


currents are sort of duals to those in (6.3) in that B E and E —> —B. The spin current 
and its relation to QED is under current study. 
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